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Abstract — We have performed holographic interferometry measurements of the dissolution of 

the (010) plane of a cleaved gypsum single crystal in pure water. These experiments have provided 

the value of the dissolution rate constant k of gypsum in water and the value of the interdiffusion 

coefficient D of its aqueous species in water. D is 1.0 x 10~ 9 m 2 s _1 , a value close to the theoretical 

value generally used in dissolution studies, k is 4 x 10 -5 mol m~ 2 s _1 . It directly characterizes the 

microscopic transfer rate at the solid-liquid interface, and is not an averaged value deduced from 

quantities measured far from the surface as in macroscopic dissolution experiments. It is found to 

be two times lower than the value obtained from macroscopic experiments. 

1 INTRODUCTION 



Pressure Solution Creep (PSC) constitutes a major plastic strain process of immersed solids under low stress. A 
crystal wetted by a solvent (generally water) and saturated by its components is at chemical equilibrium. But 
if that crystal is then submitted to an external stress, it experiences dissolution. This originates in the chemical 
potential change of the solid induced by the stress. The dissolved species then diffuse away from the high stress 
region and precipitate in a stress-free zone. This mechanism is known to play a chief role in the upper crust and 
contributes, for instance, to upper crust deformation and to the diagenesis of sedimentary rocks. 

The study of PSC in gypsum is of double int erest. First it controls th e upper crust strength in numerous 
geological situations because of its high ductility ( deMeer and Spiersl . 1997t h Additionally, PSC in gypsum also 
likely plays a role in the high ductil i ty and low strength of plaster, which is essentially made of gypsum needles, in 
the presence of moisture ( Chappuisl . 19991 ). 

To understand the mechanisms involved in PSC in gypsum, precise values of the involved quantitative parameters 
are needed. Among them, the diffusion coefficient and the dissolution rate constant are primordial. Furthermore, 
these values are also of interest in other fields. Wide areas of gypsum karsts exist world-wide The instability of 
these karsts and potential collapse undermines the importance of understanding karst e volution and subsequ ently, 
the importance of knowing the precise values involved in the water-gypsum interaction ( Jeschke et all [200lh . The 
presence of the dissolved species of gypsum in water influences the dissolution of minerals containing pollutants 
(especially i n uranium mines). H ence the study of such contamination also necessitates values regarding gypsum 
dissolution ( Kuechler et at . 20041). The existenc e of a large quantity of dissolved mineral may also alter the quality 
of the drinking water ( Raines and DewersUl997T ). Lastly, knowledge of gypsum dissolution/pre cipitation behaviour 
is req uired in the oil and gas industry where gypsum is a common annoying "scale" mineral ( Raiu and Atkinsonl 
Il990h. 



S urprisingly, no diffusion coeffi cient measurements of gypsum in water are available in the literature (jLobo and Quaresma . 

19891 : Zavtsev and Asevev . 1992 ). Recent studies of gypsum dissolution either estimate the interdiffusion coeffi- 



cient of the aqueous com ponents (jRaines and Dewersl . Il997[ ) or make use of a numerical value co mputed in 1971 



( Barton and Wilde] . 1971 ) with a Nernst-Hartley equation from the tracer diffusivity of the ions ( Jeschke et al 



20011 ). Classical dissolution measurements typically use rotating-disk setups, batch-experiments or column-experiments 
with powder or crystal samples. In these apparatuses, the concentration during dissolution is monitored by tech- 
niques such as titrimetry, conductimetry, atomic absorption spectroscopy, mass spectrometry or colorimetry. During 
these experiments, water is usually flowing and the overall measured concentration stems from a combination of 
diffusion, forced convection and dissolution. If empirical equations always enable one to describe the measurements, 
the link with pure phenomena is not straightforward and n eeds accurate knowledg e of the surface morphology 
and of the transport properties of the system to be reliable (jJeschke and Drevbrodtl . l2002al ) . To overcome these 
experimental limitations, we have carried out local, instead of global, experiments of gypsum dissolution in water 
with an alternative technique, holographic interferometry. This method allows simultaneous determination of the 
dissolution rate constant of gypsum in water at the dissolving interface and of the diffusion coefficient of its com- 
ponents in bulk water. The fundamental difference between our measurements and standard ones lies in the fact 
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that we look directly at the surface behavior instead of deducing this behavior from quantities measured far from 
the surface. 



2 EXPERIMENTS 

Historically, real-t ime holographic interferometry has been recognized as a valuable tool in solution chemistry 



([Knox et all . 119671 ). One of its advantages is that defects along the optical path (except in the working cell) com- 
pensate between the two light e xpositions of the hologram (see below), thus requiring less demanding setups than 
classical interferometry devices (jColombani and Bertl . (20061 ). Furthermore the observation of the two-dimensional 



concentration field in the liquid enables the identification of any non-diffusive fluxes (na tural convection, gravit a- 



tional instability, . . . ), thus providing a guarantee of the reliability of the measurements (|Colombani et all . 119981 ). 

Our setup is designed as follows (see Fig. [TJ. The beam of a Diode Pumped Solid State laser (A = 532 
nm wavelength) is divided into a reference beam and an object beam by a 50/50 beam splitter. The two beams 
cross a half- wave plate and a Glan- Taylor prism (vertically polarizing), which act as a polarizer/ analyzer system 
controlling the intensity of the beam. The two beams are both expanded and spatially filtered by a microscope 
objcctivc/pinholc/convcrgcnt lens set. A parallel plate installed on a rotating mount is inserted between the pinhole 
and the convergent lens of the reference beam. This device allows the removal of eventual parasitic fringes on the 
initial interferogram. The two beams interfere on the holographic plate, their polarization vectors being parallel 
and colinear to the plate. The entire setup is located on a vibration-damping structure. 

Air motion in the laboratory causes a random change with time of the path length along the beams, creating 
parasitic vibrations of the fringes. To reduce this source of uncertainty, the entire optical table is covered with 
a closed shelter and the real-time optical data acquisition is performed in a contiguous room instead of in the 
laboratory. 

The primary interest of holography lies in the fact that the interference pattern of the object and the reference 
beams contains the amplitude (as in classical photography) and the phase of the object beam. In this way, the 
amplitude and phase of the object at time to are recorded on the hologram. Therefore, when both the object 
and the hologram (containing the memory of the object at time to) arc enlightened at time t by the laser, the 
phase difference A<fi of the object between times to and t are visualized through N = A<fi/(2ir) interference fringes 
(interferograms). The interferograms are visualized and recorded with a Charge-Coupled Device camera, a monitor 
and image acquisition software. 

This device, thus, gives access to any evolution of the phase A0 = 27reAn/A (e path length in the optical cell) 
and consequently to any evolution of the refractive index An of the liquid (all other optical components being kept 
unchanged). Therefore, the change of solute molality Am = An/(dn/dm) = NX/(e(dn/dm)) in the solution can 
be registered (dn/dm is the derivative of the index of refraction n of the solution with the solute molality m). 

The crystal/solution system is located in an optical cell inside a copper structure with circulating thcrmostated 
water. The temperature of the water bath is regulated by a PID regulator, which probes the temperature with 
a platinum thermistor inserted into the copper structure. Indeed a temperature variation may be the source of a 
solubility, diffusivity and dissolution rate constant change. Therefore, constant temperature must be guaranteed 
(±0.01 K in our experiments). 

An experiment proceeds as follows: The optical cell is filled with ultrapure water. The liquid is left in the 
copper structure for a few hours in order to obtain a homogeneous temperature in the system. Subsequently, a 
reference hologram is taken. Then a cleaved gypsum single crystal (CaS04,2H20) from the Mazan mine (Vauclusc, 
France) of rough dimensions 5x5x1 mm 3 is introduced into the cell, this time being considered as the time origin 
of the reaction. Interferograms arc then recorded periodically (sec Fig. [5]). The white homogeneous zone above the 
fringes in the first interferograms is considered to be pure water (the inhomogencous grey zone is a parasitic fringe). 
Indeed, no change of color is observable in this part of the cell compared to the reference hologram, so the phase 
has remained unchanged and we may infer that no discernible gypsum has migrated yet to this zone. Proceeding 
downwards, an initial black fringe is encountered, revealing a change of tt in the phase of the transmitted beam, 
compared to the flat-colored zone just above. Hence, we can conclude that the solute molality has increased at this 
point until it has reached a value for which the resulting change of phase is ir. This corresponds to a molality change 
Am = X / (2e(dn / dm)) (cf. above). The liquid just above the fringe has been identified as pure water therefore the 
molality all along this fringe is merely to = Am. The adjacent white fringe, immediately below, corresponds to a 
7r change of the phase compared to the black fringe, and a 2-7T evolution compared to the flat-colored zone. The 
molality at this level is then m = 2Ato. The adjacent black fringe, just below, corresponds to a 7r change of the 
phase compared to the white fringe, and a 3tt evolution compared to the flat-colored zone. The molality at this 
level is then to = 3 Am. The molality in the whole cell is then gradually reconstructed (see Fig. [3]). 
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The limiting factor of the spatial resolution of the interferogram is the number of pixels in the camera sensor. 
The larger the pixel number, the smaller the area depicted by one pixel. The size of the chip of our camera is 
752x582 pixels and our resolution 40 /im. This value is small enough to provide us with a correct discretization of 
the molality in the cell. 



3 DATA ANALYSIS 

The investigated chemical reaction can be written as : 

CaS0 4 • 2H 2 0(s) ^ Ca 2+ (aq) + SO| - (aq) + 2H 2 0(1) (1) 

It takes place at the lower end of our cell. Considering the symmetry of the experiment, it can be reduced to a 
one-dimensional diffusion problem along the vertical coordinate z. In the absence of other added salts, the two 
aqueous species are present at each z value with the same molality to preserve local electroneutrality. Therefore, 
m Ca 2+ (z, t) = m s0 2- (z, t), m Ca 2+ (z, t) and m S0 2- (z, t) being the molalities of Ca 2+ (aq) and S04~(aq), respectively, 
at position z and time t. Our investigation method is solely sensitive to the change of refractive index n, regardless 
of the species causing this change. In fact, both ions induce the modification of n, plus a possible contribution of 
ion pairs and impurities. Therefore, we experimentally access an effective molality of dissolved gypsum m(z,i) m 
m Ca 2+(z,t) = m S0 2-(z,t). 

To evaluate m along the vertical coordinate as a function of time, Fick's second law is used, if the diffusion 
coefficient is considered as constant in our molality range (see Section [4]), this law reads: 

dm\ ^ / d 2 m\ 

with D as the interdiffusion coefficient of the solution, and must be solved in space and time with the following 
boundary conditions (all the symbols used are summed up in Table 1): 

• Mass balance at the dissolving interface implies the amount of gypsum leaving the solid and entering the 
liquid to being equal: 

-Pdissolution(0,£) = -Fdiffusion(0, t) (3) 

where F is the flow rate and the origin of coordinate is the geometric solid-liquid interface, considered as fixed 
(actually moving of less than one pixel in the course of a typical experiment). 

The expression of the dissolution flow rate is: 

-Fdissolution(C), t) = (4) 

with £ the advancement var iable (amount of transformed reagent). If we follow a transition states theory 
( Shiraki and Brantley . 19951 ). the rate of advancement d£/dt takes the form: 



(it 



1 ""PC" iff- 



(5) 



where k is the dissolution rate constant of gypsum in water, s r is the total dissolving gypsum-water interface, 
a so i is the activity of the reagent (conventionally taken as one for a solid), A is the chemical affinity of the 
reaction, R is the gas constant, T is the absolute temperature and is a constant. The chemical affinity of 
dissolution, in other words the change of Gibbs free energy when gypsum changes its thermodynamic state 
from solid to dissolved, is A = — RThiil. In this expression, Q stands for the supcrsaturation: 

_ a Ca 2 + Q so^- a H 2 o (6) 
K ' 

ax is the activity of the subscripted aqueous species and K sp is the solubility product of the reaction of Eq. 
The solubility product can be written as: K sp = a^ 2+ a^ 2 - (offfo) 2 : wnere a 3f* represents the activity 



'Our standard state is characterized by a unit activity for pure gypsum (solid reagent) and pure water (solvent) at any temperature 
and pressure. For the aqueous species of the solute, it is a unit activity in a 1 mol kg -1 solution representing an infinitely dilute solution, 
for any temperature and pressure. 
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of the subscripted species at chemical equilibrium. As the solubility of gypsum is small, the activity of the 
solvent au 2 o is always considered to be equal to one. If we introduce the mean ionic activity coefficient r y±, 
geometric mean of the two single ion activity coefficients, the supersaturation can be rewritten: 



n 



( 7 sat)2 m sat 2+m sat^_ 



(7) 



nix ' is the molality of the aqueous component X at chemical equilibrium, and 7^?' is the mean ionic activity 
coefficient at chemical equilibrium. Here we make the reasonable assumption that at all the ionic molalities 
achieved in our experiments, the ionic strength remains small enough so that the mean ionic activity coefficient 
does not depart from its value in pure water. This term can thus be suppressed from the expression of Eq. [JJ 

In the absence of experimental or theoretical determination of the value of the constant AT in Eq. [51 we 
have used the value successf ully used for fitting calcite dissolution measurements ( Rickard and Siobera . Il983 ; 
Shiraki and Brantley . 1995): J\f = 1/2. This choice results in the dissolution flow rate at the solid-liquid 
interface: 



dissolution 



(o,t) = ks r 1 



^m Ca 2+ (0,t)m SO 2-(0,^) , 



f Ca 2 + 



m sol~ 



(8) 



We will identify the geometric mean of the ionic molalities with our effective molality m. The geometric mean 
of the ionic solubilities is well known and will be written m sat . Finally, the dissolution flow rate at the surface 
reduces to: 

F, dissolution (0, t) = ks r (\ - • ( 9 ) 

The diffusion flow rate is linked to the diffusion flux at the dissolving interface ./diffusion (0, t) merely through: 



^diffusion (0, t) — --Fdiffusion(0,t) 
S 



(10) 



where s is the cell section perpendicular to the mass transport. This flux derives readily from Fick's first law: 

'dm(0,ty 



ffusion 



(0,t) = -Dp 



dz 



(ii) 



with p as the density of the solution. The introduction of Eq. [9lfTUland[TTlinto Eq. [3]induces for the boundary 
condition at the lower end of the cell: 



[ dz 



kj3 
~D~p 



1 



m(0,t) 



(12) 



where we have introduced P = s r /s. 



• At the upper end of the cell (water-air meniscus), the boundary condition is merely written Jdiffusk>n(£, t) = 0, 
L being the height of water in the cell. 

To get a tractable analytical expression of m(z,t), we make a further assumption concerning the geometry of 
our experiment. Our cell is high enough and the investigated time short enough to consider that gypsum diffuses 
in a semi-infinite medium. This hypothesis has been experimentally verified, paying attention to the fact that the 
highest molality fringe has no t reached the top of the cell at the end of the measurement. The solution of Equation 
[5] is in this case (jCrank , 119751) : 



m(;z, t) = 777 s 







6rfc f ) 


— exp 


\2Vmj 



kBz , kd n9 „ 
+ (TT7=) 2 ^ 



Dpm 



' Dpm s 



x crfc 



fe/3 



2VDt Dpm? 



'Dt 



(13) 



erfc is the complementary error function. The values of 777 sat , D and k are obtained through a fit of our 
experimental 777(2, t) curves with the above expression. The shift to the right of the molality curves (Fig. [3]) can 
be viewed as a signature of diffusion, and the shift of the molality curves upward can be viewed as a signature of 
dissolution. 
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4 RESULTS 



dn/dm is a value of primary importance, but unfortunately no experimental determination is available. Therefore, 
we have carried out refractive index measurements with an Abbe refractometer in pure water and water where 
gypsum was dissolved unt il saturation at ambient temperature (m sat = 15 mmol kg -1 corresponding to 2 g h 1 
(jRaiu and Atkinsonl I1990T O. These experiments were not far from the resolution limit of the apparatus. A (0.019 
± 0.007) kg mol -1 value was found, which was used for the refractive index into molality conversion. 

To guarantee the validity of the data analysis , Equation [2] requires D to remain constant with molality. Unfor- 
tunately no experimental determination of D exists, and furthermore of the dependance of D on molality. To obtain 
an initial, crude idea of the variation AD of the diffusion coefficient of gypsum in water in the molality range Am 
between pure water and saturated solut ion (Am = m sat = 15 mmol kg -1 ), we have collected t he 3D /dm values of 
other calcium and sulfate salts in water ( Lobo and QuaresmaLll989l : IZavtsev and AsevevLll992l ) and computed their 
AD for the same molality range. The worst case is ZnSC>4 at low molality with AD/D ~ 20%. Evidently, no strict 
conclusion can be inferred for gypsum from that value, but a relatively weak dispersion of D in our experimental 
cell can be expected. 

Precipitation and dissolution of minerals are known to be strongly influenced by the presence of impurities. 
To evaluate the purity of our single crystals, we carried out the Electron Probe Microanalysis of a sample of 
our material. No elements, other than the atoms constituting gypsum, have been detected with a concentration 
greater than 0.1%. The known impurities found in the gypsum of the Mazan mines, possibly present with a lower 
concentration, are dolomite rocks, quartz, anhydrite and celestite. 

For the computation of (3 = s r /s used in Equation 1131 the cross-diffusional section area s is taken as the 
horizontal section of our optical cell (s = 0.95 cm 2 ) and s r is taken as the dissolving surface area of the single 
crystal (s r « 0.3 cm 2 ). 

We have measured the pH of our solution at the end of the experiment and observe a pH of 5.7, which roughly 
corresponds to the pH of the CG^/HCO^" system stemming from the dissolution of the atmospheric CO2 in water. 

Seven measurements, all at 20.00°C, have been performed to check the reproducibility. The fit of the experi- 
mental points by the theoretical l aw of Eg . [TBI brings a solub ility m sat = 15.6 mmol kg -1 , which is very close to the 
expected value of 15 mmol kg~ : (jRaiu and Atkinsonl . [T990h . This agreement validates our analysis procedure and 
particularly the chosen dn/dm value. We find a diffusion coefficient D = (1.0 ± 0.1) x 10~ 9 m 2 s , a value very 
close to the 0.9 x 10~ 9 m 2 s _1 value derived from Nernst-Hartley equation, generally used in dissolution studies 
( Barton and Wildd . Il97ll ). We find a dissolution rate constant k = (4 ± 1) X 10 5 mol m 2 s 1 . The statistical 
uncertainty is the standard error, computed from the seven measurements. 



5 MIXED KINETICS 



There has been a strong debate to clearly determine the nature and the kinetics of the slowest step of the dissolution 
process, controlling the kinetics. If mass transport in the liquid (proportional to the diffusion coefficient D) is slow 
compared to the chemical reaction at the solid surface (proportional to the dissolution rate constant k), ions are 
rapidly unbound from the solid but slowly transported away in the solution, hence the dissolution global kinet- 
ics is controlled by molecular diffusion. Conversely, if mass transport proceeds quickly compared to the reaction 
rate, the kinetics is controlle d by the slow ion detachment from the mineral. Both macroscopic studies — rotating- 
disk, column-experi ments . . . (jRaines and Dewerd . [l997l : lDrevbrodt and Gabrovsekll200oH iDewers and Rained . 2000; 
Jeschke et al. , 2001 ) — and microscopic studies — atomic force microscopy ( Bosbach and RammenseeT 19941 : lHall and Cullen , 
19961 ) — have tried to address the question. 



One of the reasons of the diversity of the available results lies in the fact that the balance bet ween transport 



and re action strongly depends on the geometry of the system and on the thermodynamic conditions. iJeschke et al 



(|200lh conclude with mixed kinetics, where both effects are of comparable magnitude, with the linear kinetics of 
Eq. [5] (^dissolution ~ m i except very close to the saturation) for the dissolution. This conclusion is drawn from a 
combination of batch dissolution and rotating-disk experiments. Therefore, the authors make use of the formalism 
of these kinds of methods to evaluate the two phenomena: the dissolution velocity is evaluated by k s , an empirical 
dissolution rate constant, and the diffusion velocity is evaluated by k t = Z?c eq /e, a transport constant (c eq stands 
for the molarity of aqueous calcium at saturation and e for the thickness of the diffusion boundary layer surrounding 
the crystal). The authors find k t = 1.5 x 10 -3 mol m~ 2 s _1 from batch experiments and k s = 1.3 x 10 -3 mol 
m~ 2 s _1 from both experiments (for c/c oq < 0.94). These two values are so close that the mixed aspect of the 
dissolution kinetics is unambiguously assessed for this situation. Beyond the interest of this result, one can see that 
the selected evaluation quantities are strongly linked to the experimental methods. 
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Besides, IMurphv et al. ( 19891 ) have proposed the ratio of the maximum diffusive and reactive flow rates a = 
Dpm sat /((3kL) to be a dimensionless number used to discriminate between transport-controlled and reaction- 
controlled dissolution. This number has a more universal impact, with no quantity intrinsic to the chosen ex- 
periments (like e above). Pure transport control corresponds theoretically to a = 0, pure surface control to a = oo 
and mixed control to a = 1. These authors have numerically investigated realistic geological dissolution configu- 
rations of quartz and calcite in water in geometries comparable to ours (zero-flux boundary) and conclude that at 
values of a between 10 -2 and 10, the rate of evolution of the system is controlled by mixed surface reaction and 
diffusion kinetics. 

Our experimental values lead, with L = 3 cm and p = 10 3 kg m -3 , to a = 0.05. This reveals mixed kinetics, 
according to the classification of Murphy et al. At ambient temperature and pressure, the mixed kinetics of 
dissolution of gypsum in water seems to be a robust feat ure, in our geometry (a reacting surface at one boundary 
and a zero flux at the other) as well as in the geometry of Jeschke et al. ( 200ll) (a reacting surface at one boundary 
and a constant composition reservoir at the other). 



6 COMPARISON BETWEEN LOCAL AND GLOBAL DISSOLU- 
TION MEASUREMENTS 



The comparison of our dissolutio n rate constant k values with results of global measurement s requires particular 
care. As has been often stated (jRickard and Sjobergl . Il983l : Ijeschke and Drevbrodtl . l2002bl ) . derivation of rate 
constants from these kinds of experiments is a difficult task and the procedure strongly depends on the device 
geometry and sample morphology. 

To enter into the details of a global experiment, the overall flux 1Z in a vessel of volume v is computed from the 
time evolution of the overall concentration c and is considered to be identically equal to the diffusion and dissolution 
fluxes: 1Z = Jdissoiutkm = ■^diffusion = {v / s)(dc/ dt). But as water is flowing in the vessel, this equality is only valid 
at the frontier between the diffusional boundary layer and the bulk liquid, in other words along the section called s 
above. Indeed this section corresponds to the geometrical locus where mass balance of Eq. [3]applies. This statement 
remains true exclusively for mixed kinetics. For pure reaction-controlled dissolution, mass balance applies at the 
solid-liquid interface, called s r above, and for pure transpor t-controlled dissolution, no information on dissolution 
coefficients can be obtained (jJeschke and Drevbrodtl . l2002al ) . 

Between the solid surface and the extremity of the boundary layer, there is no loss of solute, therefore mass 
balance imposes the flow rates at these two surfaces to equalize. At the solid-liquid interface (of area s r ), the flow 
rate F r can be deduced from an analog to Eq. [9) F r — s r k(l — c/c cq ). At the top of the diffusional boundary layer 
(of area s), the flow rat e can be written = sk s (l — c/c cq ) where k s is the empirical dissolution rate constant 
introduced in Section [5] (jJeschke et all 120011 ). Accordingly F r = Fd implies k = sk s /s r and we are now able to 
compare our microscopic dissolution rate constant k with a rate constant k s deduced from a global dissolution 
measurement. 



Recently. I Jeschke et al. ( 2001 ) found a linear dissolution kinetics and k s = 1.3 X 10 3 mol m 2 s 1 , as previously 



mentionned. The authors estimate the geometric specific surface area of their powders by optical microscopy to 
range between 60 and 73 cm 2 g _1 , depending on the gypsum type. We consider this surface to be similar to 
our diffusional section s. The BET-surface of these powders is 1100 cm 2 g^ 1 . We identify this surface with our 
dissolving interface s r 0. Therefore, one finds a s/s r ~ 1/15 ratio, which gives a k = 9 x 10~ 5 mol m~ 2 s _1 value. 

Beside the uncertainty on the surface ratio, a factor of 2 or 3 between their global measurements (9 x 10 -5 mol 
m~ 2 s _1 ) and our local measurements (4 x 10 -5 mol m~ 2 s _1 ) seems to exist at first sight. But the large differences 
between the two methodologies could account for this discrepancy and we can consider these two values as being in 
fair agreement. One should mention that the k we have measured concerns exclusively the dissolution of the (010) 
plane, whereas the k s of global experiments is an average of the dissolution rate constants of all the reacting planes 
of the dissolving powder. Therefore, the two values are not strictly comparable. 

As a summary, microscopic measurements make the data interpretation easier for three reasons. First, there is 
no uncertainty on the surface where mass balance must be applied. Second, there is no flow in the cell and therefore 
no hydrodynamical assumptions and computations are needed. And finally the study of a single crystallographic 
interface avoid to obtain a multifaceted average dissolution rate constant. 



f Obviously this is questionable because the adsorption sites in the BET method and the dissolving sites in the dissolution experiment 
may slightly differ. But this assumption should at least bring the correct order of magnitude. 
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7 CONCLUSION 



We have used a non-invasive interferomctric method to access to the microscopic dissolution rate constant of 
a cleaved (010) surface of gypsum in water. This constant has been revealed to be two times lower than the 
same constant measured by macroscopic averaged methods, which can be considered to be in fair agreement. The 
interdiffusion coefficient of the aqueous species of gypsum in water has also been measured during these experiments 
and exhibits a value close to the theoretical value generally used for dissolution studies. Now, to deepen the 
comprehension of pressure solution creep, we plan to perform the same experiments with gypsum crystals under 
stress. The knowledge of the influence of an uniaxial pressure on the dissolution parameters may shed new light on 
the mechanisms of PSC. 
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List of Tables 

1 Description of the symbols used 



Symbol quantity units 

j3 geometrical constant (s r /s) 

7± mean ionic activity coefficient 
7j_ at mean ionic activity coefficient at equilibrium 

A<fi phase difference 

AD diffusion coefficient difference m 2 s _1 

Am molality difference mol kg -3 

An refractive index difference 

e thickness of a diffusion boundary layer m 

A laser wavelength m 

£ advancement variable mol 

p density of the solution kg m -3 

fl supersaturation 

A chemical affinity of dissolution J mol -1 

ax activity of species X 

Ox* activity of species X at equilibrium 

c molarity of aqueous calcium mol m 

c eq molarity of aqueous calcium at equilibrium mol m 

D intcrdiffusion coefficient of the solution m s -1 

e path length in the optical cell m 

F flow rate mol s _1 

J flux mol m -2 s _1 

k dissolution rate constant mol m -2 s _1 

k s empirical dissolution rate constant mol m -2 s _1 

kt transport constant (= Dc eq /e) mol m -2 s _1 

K sp solubility product 

L height of water m 

m effective molality of dissolved gypsum mol kg -1 
m sat effective molality of dissolved gypsum at equilibrium mol kg -1 

mx molality of species X mol kg -1 

rrix' molality of species X at equilibrium mol kg -1 

n index of refraction 

N number of interference fringes 

J\f numerical constant 

R gas constant J mol -1 K _1 

1Z overall flux mol m~ 2 s _1 

s area of the section perpendicular to mass transport m 2 

s r area of the dissolving interface m 2 

T absolute temperature K 

t time s 

z vertical coordinate m 

Table 1: 
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List of Figures 



Holographic setup: DPSS laser, beam splitter, Glan- Taylor prism, half-wave plate with adjustable 
polarization direction, beam expander comprising a microscope objective, a pinhole and a converging 
lens, thermostatcd optical cell, holographic plate, mirrors, parallel plate on a rotating mount, camera. 1 13 
Interferograms during the dissolution of a gypsum single crystal in pure water 9, 60, 120, 180 and 
400 minutes after the beginning of the experiment. A black bar indicates the position of the crystal 

in one of the interferograms 1 14 

Evolution of the molality with vertical position during the experiment of Fig. [21 9, 20, 30, 60, 75, 
105, 120, 180, 240 and 360 min. after the beginning of the experiment. Circles are for experimental 
points and lines for theoretical fit. The horizontal line is the solubility [15 
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